Fiche 002 Niveau : 0 


Tirage aléatoire 


1 Introduction 


L'exercice qui suit s’est présenté naturellement alors que je devais programmer un tirage équirépartie 
sur un ensemble d’entiers de la forme {0,:-- ,n — 1}. Pour cela je disposais d’un générateur pseudo- 
aléatoire cryptographique faisant un tirage au sort de manière équiprobable sur un bloc constitué d’un 
nombre exact d’octets, c’est à dire sur un ensemble de la forme {0,--: ,2% — 1}, Comment à partir 
d’un tel outil définir un tirage aléatoire équiprobable sur {0,--: ,n —1}. 

J’ai donc rédigé ce petit exercice qui donne comme solution un algorithme probabiliste (11 peut avoir 
une exécution infini, mais la probabilité d’une telle exécution est nulle et en fait il aboutit pratiquement 
très vite). 


2 Exercice 


Soient n et s deux entiers tels que n < 2°. On note P la probabilité équirépartie sur l’ensemble fini 
B = {0,---,2° — 1}. Soit a un entier appartenant à l’ensemble À = {0,--: ,n — 1} et soit : 


M ={xeBl|x modn=a}. 


1) On suppose dans cette question que l’entier n est lui-même une puissance de 2, c’est-à-dire que 
n = 2*, Montrer que tous les ensembles M, sont equiprobables. On peut donc dans ce cas se servir 
du passage au modulo pour faire un tirage équirépartie sur l’ensemble À des entiers positifs < n. 


2) On suppose maintenant que n n’est pas une puissance de 2. On peut donc trouver un entier k tel 
que : 
Der nee on 


a) Donner un exemple de deux ensembles 17, et M, qui n’ont pas la même probabilité. Donc le 
passage modulo n ne permet plus de faire un tirage équirépartie sur l’ensemble À des entiers positifs 
<n. 


b) Plus généralement montrer que pour toute fonction f de B dans À, les n sous ensembles 
f7 (a) où a est un élément de À n’ont pas tous le même nombre d’éléments. 
Ceci montre que la stratégie suivante pour tirer au sort un élément a € À : on tire au sort un élément 
b € B et on calcule a — f(b) ne conduit pas à un tirage équiréparti. En particulier, la solution de 
proportionnalité qui vient immédiatement à l’esprit et qui consisterait à tirer b au sort et à calculer un 
arrondi de 2 ne donne pas une probabilité équirépartie sur À. 


3) On propose alors pour obtenir à partir d’un tirage équiréparti sur B un tirage équiréparti sur A 
de procéder de la façon suivante : 
On tire un élément de B. S’il est < n on a terminé on a tiré un élément de À au sort. S’il est > n on 
recommence jusqu’à tirer un élément de À. 


a) Quelle est la probabilité pour que l’algorithme probabiliste ainsi défini ne se termine pas. 
a) A-t-on ainsi un tirage d’un élément de À qui soit de probabilité équirépartie ? 

b) Quelle est l’espérance du nombre de tirages à faire ? 

c) Peut-on procéder de cette façon si l’exposant s est bien plus grand que l’exposant k. 


d) À partir des résultats précédents monter une stratégie pour tirer au sort de manière équirépartie 
un élément de À, sachant qu’on sait tirer au sort de manière équiprobable un élément de b. Calculer 
l’espérance du nombre d’itération de la stratégie que vous proposez. 


3 Solution 


1) Dans cette question, n = 2 < 2%, Soit x € {0,--:,2 — 1}. En faisant la division euclidienne 
de x par 2F on peut écrire x de manière unique sous la forme : 


Gore (Open 


Si on fixe r vérifiant 0 < r < 2, le nombre d’éléments x tels que x mod 2 = y est exactement 
égal au nombre d’éléments q < tels que 0 < q 2F br < 95, c’est-à-dire aussi le nombre d’éléments 
q > O tels que 0 < q 2° < 2° (en effet si 2° < 25 alors puisque r < 2F on a aussi 2° + r < 2 x 2F — 
2FFI < 2$, Ce nombre d’éléments est 2°—* et ne dépend pas de r. En conséquence les 2° ensembles 
A, sont équiprobables. Donc la stratégie qui consiste à tirer de manière équiprobable un élément x 
dans {0,--:,2° — 1} et prendre comme résultat x mod n permet, lorsque n = 2°, d’avoir un tirage 
équiprobable sur À = {0,---,n—1}. 


2) On suppose maintenant que n n’est pas une puissance de 2. On supposera que 2*7! < n < 2. 


a) Considérons les points a = 0 et b = n — 1 de A. Alors M, est l’ensemble des multiples de 


n qui sont dans B = {0,:--,2° — 1}. C’est aussi le nombre des entiers q > 0 tels que qn < 2°. En 
conséquence le nombre des éléments de AZ, est : 
9e 
ie | | 
n 


Le nombre des éléments de M, est quant à lui, le nombre des entiers q > 0 tels que qgn + n — 1 < 2°, 
c’est-à-dire (q + 1)n < 2° + 1 où encore (q + 1)n < 2°. En conséquence : 


#M, +1 = | 
n 


Donc #M, = #M, — 1. 


b) Soit f une application de B dans À. Les images réciproques f-!({a}) où a € A forment 
une partition de B. Donc si ces images réciproques ont toutes le même nombre d’éléments { on a 
#B = 1#A. Comme #A = n n’est pas une puissance de 2, ceci n’est pas possible. Donc f ne 
peut pas être utilisée pour fournir à partir d’un tirage équiprobable d’un élément de B, un tirage 
équiprobable d’un élément de À. 


3) Tirages successifs. 


a) La probabilité de tirer un élément < n est n/2°. La probabilité de ne pas tirer un élément < n 
lors des u premiers tirages est (1 — n/2°)". Or 1 —n/2% < 1, donc la probabilité de ne jamais tirer un 
élément < n, autrement dit d’avoir une exécution infinie de l’algorithme probabiliste qu’on a mis en 
place est donc nulle. 


b) Si a; et a2 sont deux éléments de À, comme À C B et que chaque élément de B a la même 
probabilité d’être tiré, chaque élément de À a la même probabilité d’être tiré. On a donc un tirage 
équirépartie sur À. 


c) Soit X la variable aléatoire à valeurs dans {1, ---+ co} qui à une expérience fait correspondre 
le nombre de tirage qu’on a dû faire avant de tirer un nombre < n. L’espérance mathématique de X 
est du fait que la probabilité d’une exécution infinie est nulle : 


+00 
E(X) = ÿ pau, 
u=1 


où p,, est la probabilité pour que l’algorithme s’arrête après le u° tirage, c’est-à-dire : 


0 nn 
+ 1 = u—1 ; 
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(pour sommer la série qui intervient penser par exemple à la dérivée d’une série géométrique). 


d) La formule précédente nous montre que si 2° est grand devant n, 1l faudra faire en moyenne 
beaucoup de tirages 


e) Pour éviter cet écueil on peut profiter du résultat de 1) qui dit qu’on peut valablement com- 
mencer par prendre un le modulo une puissance de deux. Donc si 2%! < n < 2} on peut commencer 
par prendre le modulo 2. Ceci ramème le problème au cas où s = k. Et dans ce cas 


BONE 25 
nm 


Voici donc l’algorithme complet dans lequel on suppose 2-1 < n < 2. 
On exécute : 


1. On tire x au sort dans B. 


2. On calcule y = x mod 2f. 


jusqu’à ce que y < n. 


La moyenne du nombre de tirages nécessaires est plus petite que 2. 
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